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The present solar and atmospheric neutrino data together with the LSND results and the 
presence of hot dark matter (HDM) suggest the existence of a sterile neutrino at the eV scale. 
We have reanalysed the effect of resonant sterile neutrino conversions induced by neutrino magnetic 
moments in a type-II supernova. We analyse the implications of — v s and u B — v s {v s denotes sterile 
neutrino) conversions for the supernova shock re-heating, the detected v e signal from SN1987A and 
the r-process nucleosynthesis hypothesis. Using reasonable magnetic field profiles we determine 
the sensitivity of these three arguments to the relevant neutrino parameters, i.e. the value of the 
transition magnetic moment and the v e — u s mass difference Am 2 = ^tn 2 LSND ^ HDU 



PACS numbers: 14.60.Pq, 14.60. St, 13.15+g 



I. INTRODUCTION 



It was recognized a while ago p] that the simplest model to reconcile the existence of a hot dark matter (HDM) 
component in the Universe [Q] and possibly incorporating an explanation of the LSND results |3| with the solar Q| 
and atmospheric neutrino data || is to invoke a fourth sterile (electroweak singlet) neutrino at the eV scale. Recently 
there have been many papers on this scheme and its variants as well as phenomenology The effects of active to 
sterile neutrino transitions induced by the mixing on supernova physics have been studied in a number of papers . 

Here we also consider the relevance of such active to sterile neutrino conversion v e — > v s (v s denote sterile 
neutrino) but induced by a nonzero transition magnetic moment (TMM), similar to the transition magnetic moment 
between two active neutrino flavours |J. Such TMM leads to a neutrino spin- flavour precession (SFP) effect similar 
to the active-active case Q. The same way as the SFP effect may take place resonantly in the presence of matter 
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pPf , like MSW effect pi] , so does the active-sterile SFP. Here we adopt the general case where v e and are two 
independent Majorana particles. However, the sterile state can be a right-handed component of a Dirac neutrino. In 
this case, v s (D s ) should be regarded as a particle carrying lepton number —1 (1) in our notation. 

Although expected to be small, e.g. if the neutrino mass is introduced a la Dirac in the same way as the charged 
fermion masses in the standard electroweak theory the resulting neutrino magnetic moment is known to be very small 
fL2| |, /i„ ~ 3x 10~ 19 (m„/leV)/iB where m„ is the neutrino mass and [Ib denotes the Bohr magneton. However, 
several attempts have been made to construct various mechanisms to induce a large neutrino magnetic moment of 
order ~ 10~ n /iB p~3f . In fact this is a natural possibility in the context of radiative models of neutrino mass, such as 
that of the sterile neutrino models of refs. p]Jll,ttq]. 



Let us first recall here the previous bounds on the neutrino magnetic moments from laboratory experiments as 
well as from astrophysical considerations. Starting with laboratory limits, the upper bound on the neutrino magnetic 
moments comes from the v & e scattering experiments which give Jig ], 

< 1.8 x 1(T 1( W (1) 

This bound applies to the direct or transition magnetic moment of Dirac neutrinos, as well as to the transition 
magnetic moment of Majorana neutrinos. 

From the SN1987A neutrino observations the bounds on dipole magnetic moment of Dirac neutrinos has been 
derived by considering the helicity-flipping scattering processes such as v^e^ — > VR,e~ and v^p — > vrp inside the 
supernova core fl^l . By requiring that the i/r luminosity should be not too large in order to account for the observed 
neutrino data from SN1987A in Kamiokande II |l{| and 1MB J2fj detectors, the constraint \i v < 10~ 12 /is is obtained. 
Similar upper bound has been obtained also in a recent paper ]lq| by considering neutrino helicity flipping process 
through photon Landau damping in a dense relativistic plasma in a supernova core. It has also been discussed |T^ | 
that this bound could be improved to [i v < 10 _13 /xb by taking into account the absence of higher energy neutrino 
events which would be expected due to the spin rotation of higher energy Vr produced in the core to vl in the galactic 
magnetic field. Let us note that these discussion also applies to the transition magnetic moment of Dirac neutrinos as 
well as to the ones which connects active and sterile Majorana neutrinos if one can neglect the mass squared difference 
Am 2 in the process discussed above. However, it has been pointed out that |21| these bounds could be invalid if the 
resonant re-conversion of neutrinos take place in the supernova. 

One could also obtain limits from the arguments of excessive cooling of red giant stars due to the neutrino 
emission induced by the neutrino magnetic moment. Ref. [p2| gives the bound, 

Ih, < 3 x Kr 12 Ms, (2) 

which also applies both to Dirac magnetic moments, as well as for Majorana transition moments. 

In this paper we show that a non-zero v e - v s transition magnetic moment (TMM) would have an important 
effects in the presence of the strong magnetic fields 10 15 Gauss during a SN explosion [j] since there can be a resonant 
spin- flavour precession (RSFP) between the active and the sterile neutrino flavour, similar to what occurs for the 
active-active case. The latter has been extensively considered in the literature (2j|. We show how v e - v s resonant 
conversions induced by the active-sterile TMM enhance the r-process nucleosynthesis in all of the mass range relevant 
for the neutrino hot dark matter scenario || . On the other hand we estimate the restrictions on neutrino active-sterile 
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TMM from SN1987A data, and investigate the influence of these v e - v s TMM-induced conversions in the shock revival 



problem. For previous related work see refs. 17 2l|j2qj27| . In particular, ref. |27j has considered v e - v s TMM-induced 
RSFP in the context of SN physics for the case of random magnetic fields, in which case the conversion occurs in 
a non-periodic regime. Here we are considering the case of regular fields, and we discuss novel issues related to SN 
physics which were not considered in [ ^7f . Throughout this paper we assume \i v < 10~ 12 /jb so that the neutrino 
helicity flipping process in the core |l7],[l8|], discussed above can be neglected (see the discussion in Sec. IV for the 
case if \x v > 1CT 12 /z_b). 

In Sec. II we briefly describe the picture of the neutrino propagation in matter and the resonant v e — v s 
conversion JlT| . Sec. III. A discusses the implications of v e — v s conversions for the neutrino re-heating mechanism. 
In Sec. III.B we analyse the impact of our scenario in the later epoch of supernova evolution (few seconds after the 
core bounce) for SN (anti-) neutrino detection rates (Sec. III.B), as well as for r-process nucleosynthesis (Sec. III.C). 
In Sec. IV we summarize our results. 



II. THE ACTIVE-STERILE NEUTRINO RESONANT SPIN PRECESSION 

In our discussion we only consider the v e — > v s and v e — > v s conversion channels, where v 8 (y s ) is a sterile 
neutrino, due to resonant spin precession. For simplicity we neglect in what follows the mixing between v e and v s so 
that the 4x4 evolution Hamiltonian for the neutrino system || reduces to an effective 2x2 system. Moreover we 
consider the Am 2 = m 2 s — m 2 > case, i.e., sterile state is heavier. The evolution of the v e — v s system in the matter 
background with non-zero magnetic field is determined by the following Schroedinger-like equation, 



dr \ v s ) \ [i v B J \ v, 

where [i v is the neutrino magnetic moment and B is the magnetic field strength perpendicular to the neutrino 
trajectory. The effective potential V e for v e arises from the coherent forward neutrino scattering off-matter constituents 
p| and is given by 

V. = ^(F„-ir w ) = ^(3n-l), (4) 
rriN 2 ZrriN 

n e 

Y e — ■ , Y n 1 Y e . 

n e + n n 

Note that there is no potential for v s , i.e., V s — 0. Here Gf is the Fermi constant, p is the matter density, mjv is the 
nucleon mass and n e and n n are the net electron and the neutron number densities in matter, respectively. Note that 
charge neutrality n p — n e is assumed. For the v e — > v s system the matter potential just change its sign. 



The resonance condition is: 



Let us note that for Am 2 > 0, either v e — > v s (for V e > i.e. Y e > 1/3) or v e — > v s (for V e < i.e. Y e < 1/3) 
conversions take place. This is important because, as we discussed in ref. [fjj, in the region above the neutrinosphere 
the matter potential V e changes its sign due to the different chemical content. 



2 It is reasonable to neglect the contributions from the neutrino background to the effective potential, since the neutrino 
densities in the relevant regions are small. 
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FIG. 1. Typical \ V e \ profile versus the radial distance r at t < 1 s after the core bounce. The solid and dashed lines correspond 
to positive and negative potential, respectively, and the position where V e = is denoted by r* . The Am 2 values for which a 
E — 10 MeV neutrino undergoes resonant conversion, for the \V e \ value on the left ordinate is indicated in the right ordinate. 
This is relevant for the shock revival considerations. 
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FIG. 2. Same as Fig. 1 for t > 1 s. This is relevant for the r-process and v e signal considerations. 



In Figs. 1 (earlier epoch) and 2 (later epoch) we reproduce the typical profiles of corresponding \V e \ in the 
region of interest, taken from rcf. J^]. They are obtained from the Y e and p profiles taken from Wilson's supernova 
model outside the neutrinosphere, at two different times, t < 1 s and t > 1 s post-bounce. There is a point where Y e 
takes the value 1/3 (i.e. V e = 0), as indicated by r*. This position corresponds to r* as 160 km and 12 km for the 
earlier and later epochs, respectively. Clearly, the effective potential V e changes its sign from negative to positive at 
the point r*. 

The resonance condition in eq. (|5|) provides the Am 2 value for which neutrinos with some given energy can 
experience the resonance for a certain value of the potential (or equivalently, at some position r) . For convenience, in 
the right ordinate of the Figs. 1 and 2 we have also indicated the corresponding Am 2 values for which neutrinos under 
go resonance, for typical neutrino energy E = 10 MeV. We see that for Am 2 > 10 2 eV 2 only v e — > v s conversions can 
occur, and this happens in the region where Y e < 1/3. On the other hand for Am 2 < 10 2 eV 2 three resonances may 
occur. The i? e 's are first converted, say at r± < r* , then there are two resonance points at r-i (~ r\) and r^, where 
^3 > f2 > f* (i-e. in the region where Y e > 1/3), for the v e — * v s channel. A schematic level crossing diagram for 
v e — v s and v e — v s system may be seen explicitly in Fig. 3 of ref. jq] . Note that the v s 's originated from the first v e 
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conversion (at r 2 ) can be re-converted into v e 's at the second resonance (at r 3 ). 

In what follows we will employ the simple Landau-Zener approximation p8| . ^9| in order to estimate the survival 
probability after the neutrinos cross the resonance. Under this approximation, the v e (or P e ) survival probability is 
given by 



P = exp 



exp 



2 LI 



-1.1 x 10" 



dV e km 
dr eV 



(6) 



where L™ s is the neutrino oscillation length at resonance. Notice that for Ar/L™ s > 1 the resonant neutrino conversion 
will be adiabatic pJ] |. 

As we discussed in ref. ||, depending on the Am 2 values, the v e — v s system may encounter two resonances. 
In order to take this into account we compute the v e survival probability after the second resonance as follows, 



P(Ve 



s ) = P(r 2 )P(r 3 ) + [1 - P(r 2 )][l - P(r 3 )], 



(7) 



where P(r2) and P(r 3 ) are the survival probabilities calculated according to the eq. 
resonance positions r 2 and r 3 , respectively. 



at the first and second 



In this work we will neglect the spin precession due to the galactic magnetic field discussed in ref. JL7J! since 
for the relevant parameters of interest to us here the condition [i v Bq <C Am 2 /E, where Bq ~ 10 -6 G is the galactic 
magnetic field, is always satisfied so that the spin precession in the galactic field is strongly suppressed. 



III. IMPLICATIONS FOR SUPERNOVA PHYSICS 



Following ref. [g| we now consider the possible impact of active-sterile neutrino conversions on different aspects 
of supernova astrophysics. We analyse processes taking place at the early epochs as well as at the cooling stage, in 
order to get a feeling for their sensitivity to the underlying neutrino mass-square difference, magnetic moments and 
magnetic fields. In the following we consider only epochs after the core bounce and the neutrino evolution in regions 
outside the neutrino-spheres, where resonant conversions take place. 



A. Shock Revival 



Here we are considering the earlier epoch t < 1 second after core bounce. In the delayed explosion scenario 
pO| , ^H the neutrino energy deposition, occurring between the neutrino sphere and the site where the shock is stalled 
(about 500 km or so from the core) can re-start the shock and power the explosion. The disappearance of either v e 
or v e due to the resonant spin precession into sterile states would reduce the energy deposition rate. The neutrino 
energy deposition rate R at the stalled shock is defined as follows, 

R = KKM + YyKv^) (8) 

Y n E UBn (t)+Y p E QeP (t) ' 

<P°(E v )dE v , E' vN (t)~ J E v tr vN P{E v )<lP{E v )&E v . (9) 
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where the primed quantities Y£ and stand for the proton and neutron abundances calculated in the presence 
of active-sterile neutrino conversions. The un-primed ones correspond to the case where no conversion occurs. As 
an approximation, we have assumed, in eq. (Q) that the neutrino energy spectra cf>°(E u ) are Fermi-Dirac with zero 
chemical potential but with different characteristic temperatures T v which depend on the neutrino flavour. As in our 
previous analysis we choose T„ for v e and v e such that the typically predicted average energies, (E Uc ) = 11 MeV and 
(Ep c ) = 16 MeV are obtained in the absence of any conversion. We assume a magnetic field profile B(r) = Bo(ro/r) n 
where Bq = fO 16 G, n = 2, ro = 100 km and r is the distance from the center of the star. We have also studied the 
case where n = 3 and obtained the similar results. 
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FIG. 3. Contour plot of the reheating rate R in the Am 2 — /i£?o plane. The rate is normalized to 1 in the absence of any 
kind of neutrino conversion. Here /j,b denotes the Bohr magneton. 

In Fig. H we plot the iso-contour of R in terms of (Am 2 , ^Bq). We show the values of /i„-Bo in units of \ib • 10 16 
Gauss, where 10 16 G (at r = 100 km) might be the maximally conceivable value |32|. This is taken as a reference 
value in order to illustrate the maximal sensitivity to the neutrino magnetic moment (though such a value of magnetic 
field may affect the supernova dynamics and self-consistent calculations would be required for its justification p3fl ). 
Thus we can say that if the neutrino re-heating is essential for successful supernova explosion the parameter region 
right to the curve, say R = 0.5, is disfavoured. 

B. Implications for the detection of SN1987A u e signal 

We arc now in the later epoch t > 1 second after core bounce, in the so-called Kelvin-Hclmholtz cooling phase. 
The observed v e events from the supernova SN1987A in Kamiokande II and 1MB [^| detectors, 11 and 8 events, 
respectively, are in good agreement with the theoretical expectations. Therefore, any significant conversion of D e 's 
into a sterile neutrino species would be in conflict with this evidence. We consider the effect of active-sterile neutrino 
conversions both on the v e signal in order to analyse the possible restrictions on neutrino parameters, the neutrino 
magnetic moment and the mass-squared difference. 

We plot in Fig. ^ three contours of the v e survival probability (P) , which is averaged over the Fermi-Dirac 
energy distribution, for the v e — > v s conversion, in the (Am 2 , fiBo) plane. In the plot the left, the middle and the right 
lines correspond to (P)=0.9, 0.5 and 0.1, respectively. If wc assume that the successful observation of the SN1987A 
signal implies that at least 50 % of the expected D e signal has been detected, one can conclude that all the portion 
right to the curve P = 0.5 is ruled out. We note that due to the fact that the potential is much steeper in the 




G 



later epoch than in the earlier epoch, the bound we obtained in Fig. 4 is much weaker than the sensitivity we have 
displayed in Fig. 3. 




logQiB ( /ji B 10 ,0 G) 

FIG. 4. Contour plots of the survival probability P (figures at the curve) for the Ve. — > v s conversion at t > 1 s after the 
bounce. The region to the right of the curve corresponding to P = 0.5 can be excluded by the observation of the SN1987A i/ e 
signal. We have assumed the magnetic profile to be B(r) — B (r /r) n where ro = 10 km and n — 2. 



C. Implications for r-process Nucleosynthesis 



Here, we consider the implications of resonant neutrino conversions for the supernova nucleosynthesis of heavy 
elements [pip] . It has been proposed that the supernova could be the most promising site for heavy elements (mass 
number A > 70) nucleosynthesis, the so-called r-process p5| , in which neutrinos play very important role (36] . A 
necessary condition required for the r-process is Y e < 0.5 in the nucleosynthesis region. The Y e value in the r-process 
site is mainly determined by the following neutrino absorption reactions, 

v e + n — > p + e' (10) 
9 e +p^n + e + . (11) 

Therefore, in the nucleosynthesis region we can write Y e as follows |j4|, 

Ye ^ i + Ki^JJp^ e (E Vc y (12) 



From this expression, we notice that the v e — > v s conversion leads to an increase of Y e , whereas the v e — > v s conversion 
acts in the opposite way. As we described in detail in ref. J8|, depending on the Am 2 range, one channel dominates 
over the other one, and Y e can be increased or decreased. Using again the Fermi-Dirac energy distribution to 
average the neutrino absorption rates, we have calculated the electron abundance Y e at the site where heavy elements 
nucleosynthesis is expected to take place as functions of (Am 2 , /i£?o)- 

We present our result in Fig. 5. For a successful r-process, the region above Y e > 0.5 is ruled out. On the other 
hand we find that the supernova nucleosynthesis could be enhanced in the region enclosed by the contour Y e = 0.4, 
similar to the results obtained in ref. This region is delimited by Am 2 < 10 2 eV 2 and 10~ 16 < /i„ < 10~ 15 in units 
of [Ib for Bq = 10 16 G. Inside the contour of Y e — 0.33, Y e might get stabilised to 1/3 due to some feedback effect 
which we discussed in detail in ref. [|| . We see from the plot that the most promising mass range for the neutrino hot 
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dark matter scenario 0, Am 2 < 10 eV 2 , lies inside that where the r-process nucleosynthesis is enhanced. Moreover, 
it is neither in conflict with the re-heating process (see Fig. 3), for ix v Bq < (10 _15 /xb) • (10 16 G), nor with SN1987A 
observations (see Fig. 4). 

c~ 4.5 

> 

c\j 

| 3.5 
o 3 
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1.5 
1 
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-17 -16 -15 -14 

log((xB /|^ B 10 16 G) 

FIG. 5. Contour plots for the electron fraction Y e (figures at the curves) taking into account v e — > v B and v E — > v B conversions 
at t > 1 s after the bounce. The region to the right of the solid line labelled 0.5 is ruled out by the condition Y e < 0.5 necessary 
for r-process nucleosynthesis to occur. For the parameter region inside the Y e = 0.4 dotted contour r-process nucleosynthesis 
can be enhanced. 




IV. CONCLUSIONS AND DISCUSSIONS 



In this paper we have considered the effect on supernova physics of having RSFP conversions induced by 
non-zero active-sterile TMM. We have investigated the sensitivity of the shock revival and r-process arguments, as 
well as estimated the restrictions from the observed SN1987A D e detection rates. We have analysed the effect of 
RSFP conversions involving v e to v s or D e to v s in the region above the hot proto-neutron star in type II supernovae, 
assuming the simple magnetic field profile, B(r) = Bo(ro/r) n (n = 2,3) where Vq is taken to be 100 and 10 km for 
earlier and later epoch, respectively. If the mass of the v s is in the cosmologically interesting range and the product 
of the TMM by the magnetic field just above the neutrinosphere, /i„Bo, is larger than ~ (10 _15 /iB) x (10 16 G), then a 
significant fraction of v e and v e would be converted into v a and D s , respectively, in the region outside neutrinosphere. 
Such conversion could lead to the depletion of v e and v e fluxes, resulting in a suppression of the expected v e signal 
in underground terrestrial detectors, in contradiction with the successful observation of the SN1987A v e signal in the 
1MB and Kamiokande detectors. Hence, on this basis we can constrain the neutrino mass and neutrino magnetic 
moment by requiring that the total v e flux during the thermal neutrino emission epoch should not be significantly 
depleted by v e — > v s conversion, and can rule out the range fi^Bg > (10~ 15 /is) x (10 16 G). 

The TMM-induced v e - v s RSFP would also suppress the neutrino re-heating behind the stalled shock. We 
have found that for /i„i?o ^ (10 _16 ^s) x (10 16 G), the energy deposition by v e and v e absorption reactions during the 
shock re-heating epoch (t < 1 s after the bounce) could be significantly decreased. From the r-process argument, for 
the parameter range Am > 100 eV 2 , where the Ve. — > v s conversion is dominant, Y e at the nucleosynthesis site could 
become larger than 0.5 and hence r-process would be forbidden, leading to a disfavoured range, fi^Bo > (10~ 15 /is) x 
(10 16 G). On the other hand, for Am 2 < 100 eV 2 , r-process nucleosynthesis could be enhanced due to the decrease 
of Y e down to the minimum value 1/3 due to the efficient v e — ► v s conversion if the parameters are in the region 
(10- 16 /^) x (10 16 G) < ti u B < (10- lf >s) x (10 16 G). 
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Throughout the above discussion we have neglected, for simplicity, the production of the sterile state due to the 
neutrino helicity- flipping processes in the core, which was discussed in refs. This requires \i v < 10~ 12 /j,b- Let 

us finally discuss briefly the case /Lt„ is larger than 10 _12 /xb (but smaller than the laboratory limit (1)). As discussed 
m refs. @Hl, in this case the sterile states can be copiously produced by the helicity-flipping scattering processes in 
the SN core and they can escape freely to the outer region with higher energy. 

For the earlier epoch, relevant for shock reheating, if the parameters (Am 2 , /iBq) are in the region where the 
RSFP is efficient, the reheating rate would be increased pl| due to the re-conversion of v s — > v e and/or v s — > v & since 
the v e and v e would have higher energy compared to the standard case. Therefore, we would not get any disfavoured 
region but a positive effect from this conversion. However, for the discussion on the v e signal we still would get 
constraints due to a different reason |37j . In this case we would not have the reduction of the v e signal but instead 
we would get an increase in the expected number of events at underground detectors. This enhanced v e signal would 
follow from the increase in the average energy due to the v s — > v e re-conversion. This would allow us to disfavour a 
similar range of parameters as in Fig. 4. Finally from the r-process argument, we would have a different conclusion 
from what we obtained in Sec. Ill C. Roughly speaking, the allowed region (Y e < 0.5) in Fig. 5 would now become 
disfavoured region whereas the disfavoured region (Y e > 0.5) in the same plot would now become the allowed one. 
The reason is the following. For example, if the parameters are in the region where Y e < 0.4 in Fig. 5, only the 
resonant conversion at r% (defined in Sec. II) is adiabatic due to the fast variation at r-i but the slower variation of 
the potential at (see Fig. 2). In the anti- neutrino case, the v s — > v e conversion at T\ can be neglected, as one can 
see from Fig. 4, for magnetic moments smaller than 10 -15 hb or so. This implies that, instead of the net effect of 
the v e — > v s conversion expected in the case we considered in Sec. III.C, now we expect efficient v s — > v e conversion, 
which makes the average energy of the resulting v e higher than in the no-conversion case. This has the result of 
driving the Y e value larger than 0.5 (see eq. (12)). On the other hand one can expect the opposite behaviour for the 
region where Y e > 0.5 (and larger magnetic moment) the conversions at all the points r, are expected to be adiabatic 
and therefore important. Thus, roughly speaking, the net effect is the reduction of v e flux, when neutrinos reach the 
position where the r-process is occurring. The net effect is that the v e flux is decreased and the v e is increased. This 
way we expect that the value of Y e will become smaller than in the no-conversion case. Thus (assuming that v s and 
v s are copiously produced in the core) the region with Y e > 0.5 would now become allowed by the r-process criterion. 
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